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1. INTRODUCTION 

In the classic stable layer case in which density is a linear tracer of the mixture composition between two 
layers of fluid, turbulent entrainment destroys the available turbulent kinetic energy or it re-stratifies the flow during 
the decay of a turbulent event. This is because the stratification inhibits the vertical motions of the interface which 
marks the edge of the turbulence. In contrast to the classic case, the evaporation of the cloud droplets for which 
unsaturated air is entrained into a cloud causes cooling, and hence an increase in density, such that this additional 
downward movement due to buoyancy reversal may affect the dynamics of entrainment and instabilities of cloudtop. 

The primary interest in this study is to investigate the effect of buoyancy reversal and Richardson number 
on the simulation of turbulent entrainment into marine stratocumulus clouds or the Weddell Sea. The key question 
to be addressed in this note is: What is the effect of buoyancy reversal (evaporative cooling) on entrainment rate ? 

The oscillating grid-turbulence experiments, in which the turbulent kinetic energy is mechanically produced 
at a surface region on scales much less than the depth of the layer, have been used to investigate the entrainment 
mechanisms across an inversion for a couple of decades. For instance, many experiments have been made at 
observing and parametering the entrainment across a density stratification after the basic laboratory grid-oscillation 
experiment of Rouse and Dodu (1955). Turner (1968) showed that the turbulent motions producing the entrainment 
were directed toward the interface and the entrainment rate was a function of Richardson number in linearly stably 
stratified fluid. Oscillating grid-turbulence measurements (Thompson & Turner, 1975; McDougall, 1979; Long & 
Fernando, 1983; E & Hopfinger, 1986) indicated a spacial decay of the horizontal root mean square (r.m.s.) turbulent 
velocity w near the interface of the form waz'l where z was measured from a virtual origin. 

In order to make the comparison of entrainment rate with and without buoyancy reversal, a series of grid- 
oscillation experiments with linear density change of the mixtures (water - salt water system) was first completed 
using the same flow visulization techniques as before (Part I). In this note we organize our results using three 
parameters, namely, buoyancy reversal parameter (D), Richardson number (Ri) and Reynolds number (Re). Only one 
variable can be changed at a time. Conditions may be controlled. Most importantly, runs may be repeated. 

2. EXPERIMENTAL METHODS 


2-7 The Apparatus 

The apparatus is sketched in figure 1. A transparent tank of 28x28x60 cm height is separated into two 
compartments by a thin horizontal sliding stainless plate of 0.07cm thickness. Before a run, the compartment below 
the plate was filled with water and that above the plate with a mixture of alcohol and ethylene glycol for nonlinear 
case, using the same fluid sytem (Part I: Fig.l). It consists of a grid of 1.12 cm square bars, aligned in a square array 
of mesh size M=5.60cm. The grid, placed horizontally in the tank, can be oscillated vertically with varied frequencies 
from 1 to 10 Hz by a speed-controlled motor. The stroke was fixed at S=2cm for all experiments. The grid was 
mounted on a 0.64cm diameter connecting rod which was supported by a 7.62cm long linear bearing and connected 
with the motor by two ball bearings on a cam. Compared to the earlier non-buoyancy-reversing experiments in 
which the density of mixtures is a linear function of the concentration, the improvements here are the following. 
First, there is no interference between the connected rod and the interface during the grid-oscillation, because the rod 
dose not penetrate the interface. Second, the sliding plate sealed by o-rings is provided across the center so that the 
tank can be filled conveniently while maintaining the purity of the two layers before the run. Without this feature, 
the nonlinear density changes can cause considerable mixing during the filling even before the mechanical stirring is 
begun. Therefore, a thin removable partition between two layers is necessary for the buoyancy reversal case. An 
electrical timer recorded the elapsed time for these experiments. 

2-2 The Characteristics of Oscillating Grid Turbulence 
According to E & Hopfinger (1986), the empirical relations for calculating the r.m.s. horizontal component 
of turbulent velocity w and the integral length scale d of the turbulence near the interface are 

co = CiS 3 ^M 1/2 fz -1 and 5 = C 2 Z (1) 
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where Cj » 0.3 and C 2 = 0.1, using the same constants as Mory & Hopfinger (1985); M , S and f are the mesh 
size, stroke and frequency of the oscillating grid respectively. To reduce the wall effect, the grid was designed by 
assuming the walls were planes of symmetry as suggested by E & Hopfinger (1986). Also, the maximum frequency 
used in this study was 6 Hz to avoid the unwanted circulating motions at high frequencies (McDougall, 1979). 


(A) Experiments for linearly mixing case 

Figure 2 is the schematic diagram of the experiments in the linear case. Choosing the Z* axis as positive- 
upward (figure 2a), the density p(Z*) is given by 

p , (P> > z > (2) 

p (-Ho < Z* < z -h) * 

where h is the interface thickness, which is generally much less than z. Here 

P « Pl ^go + pogorM o 

1 -H 0 z-H 0 

and (3) 

Ap - p -pi * (po-pi)*^o±f!Q , 
z + Ho 

where Z<) is the initial position of the interface measured from the mid-plane of the oscillating grid and p 1 and po 
are the initial densities of the upper and lower layers respectively. 

For these linearly stratified experiments, the initial normalized density difference between two layers ((pQ - 
Pl)/Po) was typically set around 0.5% to 3%. The Reynolds number based on the grid’s mesh size and frequency can 
be up tol()4 and the typical Reynolds number wd/n near the interface was = 120, which may be large enough for 
viscous effects to be negligible. The mixed-layer depth z , as shown in figure 3a, is a function of time. The 
evolution of growing interface was recorded to establish a depth-time relation. The entrainment rate and Richardson 
number are then denoted as follows: 


“ <> “ 2 (4) 

where w e = dz/dt and w , d , p and Ap are given in Eq. (1) and (3). It is important to note that the scales used 
to define Rj are the r.m.s. horizontal component of turbulent velocity w and the integral length scale d of the 
turbulence near the interface instead of the velocity of oscillating grid and layer depth, using the same empirical 
relations found by E & Hopfinger (1986). 


(B) Experiments for nonlinearly mixing case 

Figure 3 shows the schematic diagram of the experiments in nonlinear case. Three different grid locations 
with respect to the interface were used: above (A), below (B) and center (C) on the interface. For these experiments, 
the stroke was set at 2 cm with a frequency of 1 - 6 Hz. The density p(Z*) as showing in figure 3b is given by 

pCZ*) - f Z * >z (5) 

where H’ is the actual mixing region indicated by color. H’ is equal to H 0 after the full lower layer has turned 
dark red. The density jump across the inversion is then approximately given by 
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- before saturation 
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1-p* z +■ H„ (6) 

where p* is the maximum density at mixing fraction p* on the buoyancy reversal curve as shown in figure lb (Part 
I). The buoyancy reversal parameter D is then denoted by 

P‘*P 1 r 

^ — — - - - - before saturation 

0 after saturation 


(7) 


We used the same empirical equation (Eq. 1) to determine the characteristic turbulent velocity and length 
scales near the inversion for the nonlinear case because the same grid-turbulence was used as that of the linear case. 
From Eq. (4) and (6), the entrainment rate and the Richardson number were then calculated. A pH sensitive dye 
identified regions of mixing was used to detect optically how 'buoyancy reversal' created by mixing process disrupted 
the inversion leading to plume-like runaway entrainment. 


3. RESULTS FOR ENTRAINMENT EXPERIMENT 


5-7 Stably Stratified Flow Without Buoyancy Reversal 
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Figure 4 shows the evolution of the flow which developed in experiment (A) at different times. At time = 7 
see after the motor started (Fig.4(a)), the stratification flattens the eddies generated from the oscillating grid. The 
eddies recoiled back, entrained and mixed fluid, turning it dark red into the turbulent flow, where Ri is about 47. As 
time increases, the mixed-layer deepens and the Ri near the interface increases, corresponding to a low value of 
entrainment rate. The experimental data for the mixed-layer depth v.s. time can be best fiued by a power law relation 
of the form z = at*\ with a being a constant for a given oscillating frequency and b « 0.18 for these experiments. 
Both constants are determined by a small least square program. Compared to the earlier linear experiments, the 
exponential constant b is close to b * 1/5 (E & Hopfinger, 1986) and b « 2/11 (Fernando & Long, 1983). Figure 5 
shows the logarithmic plots of the entrainment rate E as a function of local Richardson number Ri. For the different 
grid positions, these experimental data approximately obey 

E = C R;' 3 / 2 . for Ri>6, (8) 

where C « 2.0. For comparison, E & Hopfinger’s data are included in Fig.5, where C ~ 3.8 (E & Hopfinger, 1986), 
and C « 2.5 (Turner, 1973) with the same slope for Ri>7. It is important to note that the entrainment rate tends to 
depend less strongly on Ri and flatten to a somewhat constant level at small value of Ri(<6). The transition of Ri 6 
is close to the value of maximum internal wave generation (Carruthers & Hunt, 1986) and the E & Hopfinger’s data. 
For Ri _ 1, the vortex or eddy has the excess kinetic energy such that it engulfs fluid into itself as it rotates. As Ri 
increases (Ri>6), the buoyancy becomes important, so that the entrainment is dominated by an eddy-rebounding 
process (Linden, 1973). 

3-2 Stratified Flow with Buoyancy Reversal 

(A) Experiment A 

(1) Di=1.0 

Figure 6 shows the evolution of the nonlinear experiment A for Di=1.0 at four different times. The 
mixtures are either trapped by the interaction among these eddies or created by internal wave breaking as interface 
tilted, where upon they descend and detach slowly from the inversion, while the turbulent layer remains pure and 
unmixed (colorless). At time * 5sec after the grid started, the heavy mixtures descended, while the interface remained 
flat where Rj = 97 and D = 0.9. As time increases, more entrained fluid sinks into the non-turbulent layer, bringing 
the interface even closer to the oscillating grid, where the Ri and D are reduced. The changes in interface’s structure 
with varying Richardson number are also shown in figure 6. At about 20 minutes after starting grid the upper 
turbulent layer starts to turn red because of entraining fluid from below, where the growing interface is about 5 cm 
away from the middle plane of oscillating grid and Ri = 6. Some mixtures still continue to descend into non- 
turbulent layers due to buoyancy reversal. The key point regarding the interfacial instability is that the system is 
stable in the sense of no runaway entrainment occurrs in this case. For Dj<1.0, we saw that the turbulent kinetic 
energy of the mixed parcel created by the mixing-induced buoyancy reversal was not large enough to ruin the 
interface leading to positive feedback entrainment. Indeed, the heavy mixture only re-distributes itself into the lower 
layer (the simulated clouds) without much perturbance of the interface. The entrainment was modest during a run. 

(2) Di=2.5 

Figure 7 shows the evolution of experiment A for Dj=2.5 at five subsequent times. At t * Msec after 
starling the grid motion, the interface remained stable even at D * 1.9, when Ri *15. At t « 33sec, there was more 
dilution of the lower fluid layer by the upper fluid, reducing both the Ri and D. At t = 45sec, the interface became 
more convoluted such that runaway entrainment was incipient, while the upper layer (grid-turbulence side) still 
remained pure and unmixed, when Ri * 6 and D * 1.3 as shown in Figure 7c. A strikingly vigorous entrainment then 
occurred that disrupted the interface. Very quickly, the two layers were mixed together. The instability is apparently 
parameterized by the Richardson number as well as the buoyancy reversal parameter. 

(B) Experiment B (Di=2.0, f=4Hz) 

Figure 8 shows the evolution of the flow which developed in experiment B for an initial value of D=2.0 at 
four different times. At t = 5 sec, the interface became strongly convoluted. The associated rapid entrainment 
occurred. This rapid entrainment was sustained only a couple seconds. The heavy mixtures, upon descending, were 
quickly homogenized by grid turbulence which in turn restrained the vertical motions of the interface because of the 
strong stratification. At t =13sec, the heavy mixed fluid reached the bottom of the tank, while the interface relaxed 
back to horizontal, associated with a low entrainment rate. At the turning point in which rapid entrainment relaxed 
to slow entrainment, the measured Ri was again about 6 and D>1. As time passed, there was more dilution of the 
lower fluid layer (the simulated clouds) by the upper fluid (the simulated dry air), decreasing D and increasing Rj. The 
variations in the structure of the interface are shown in figure 8(a-d), which are photographs taken at Ri * 5.0(a), 
8.8(b), 16.1(c) and 30.3(d). 

3-3 Instability Condition 

In figure 9 we have plotted the buoyancy reversal parameter D as a function of local Ri for two series of 
experiments A and B. These experiments reveal three main points. First, it is surprising that the system is stable for 
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strong perturbations if buoyancy reversal is comparable to or less than the initial stratification (D<1.0) as shown in 
figure 9. Second, the interface can be stable for D>1.0, if Ri is much larger than 6. Third, the system is unstable in 
the sense of the inversion becoming largely convoluted associating a rapid positive feedback entrainment, if D>1 and 
Ri<6. 

3-4 Entrainment rate as a function of Ri and D 

For simplicity, in figure 10 we only include two runs of experiment (B) with two different initial values of 
D, which indicates the nonlinear entrainment rate as a function of Ri and D. It shows that the entrainment rate in the 
nonlinear case only increases about twice that of the linear case due to buoyancy reversal at D = 0.5 and Ri « 10. As 
Ri increases, the stronger stratification inhibits the vertical motion of the interface, followed by inhibiting the effect 
of buoyancy reversal. The current experimental data suggested that the effect of buoyancy reversal on the entrainment 
rate can be negligible, if Ri is large (Ri > 50) and D is small (D < 0.5). 

4 . DISCUSSION AND CONCLUSION 


A stratified interface is stable to the buoyancy reversal instability for surprisingly large values of D. A new 
instability mechanism is proposed, which considers the mixing process at the interface. For the type of density 
curves studied here, under strong perturbations, the mixed parcel must have a buoyancy reversal comparable to the 
initial stratification before the interface is unstable. This is in accord with a simple model of the interface mixing 
process, as well as aircraft observations of long-live marine stratocumulus clouds. These clouds' remarkable 
longevity in the face of finite D indicates that they can be stable (Hanson, 1984; Albrecht et al„ 1985; Siems et al„ 
1989). The present work suggests that buoyancy reversal as well as the disturbance must be large for Cloudtop 
Entrainment Instability. The effect of buoyancy reversal (evaporative cooling) does not always enhance the 
entrainment rate over that in the inert case, but it may be negligible if Ri is large (Ri > 50) and D is small (D < 
0.5). This work may shed some light on the fundamental mechanism of the breakup process of the subtropical 
stratocumulus clouds into tradewind cumulus. These results may also be related to the instability in the Weddell Sea 
off of Antarctica. 
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Figure 2: Schema tic diagram of the experiments in linear case. Three 
different grid positions with respect to the interface were used: (a) 

Below (B); (b) Center (C); (c) Above (A) the interface. 


P(Z*) 


Figure l:Sketch of the apparatus. A: Mixures of alcohol and ethylene 
glycol. B: water. It consists of a square array grid, driven by a speed- 
controlled motor. 
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Figure 3:Schematic diagram of the experiments in nonlinear case: (a) 
Experiment (A); (b) Experiment (B); (c) Experiment (C). 


238 


l 





Figure 7:The evolution of experiment A for D, = 2.5. Where Zq - \Q 
cm. f s 4Hz uid S = 2 cm: (a) I = 14.1 sec, = 15.2, D = 1.9 (b) i = 
33.6 »ec, Rj = 10.3, D = 1.7 (c) t = 45.4 iec. R, = 6.3. D = 1.3 (d) i = 
50.8 iec (e) t = 55.6 iec. 


Figure 4:The evolution of the expenmcnti for linearly mixing caie. 


Experiment A : Zq = 10on, S = 2an and f = 3Hz. (i) t = ll.liec, Rj 

* 47 (ii) t = l,840sec, Rj m 120; (b) Experiment B : Zq - 6cm, S = 

2cm and f = 3Hz, (i) t = 8.5sec, Rj - 15 (ii) t = 900.3sec, Rj - 35; 
(c) Experiment C : Zq = 0cm, S = 2cm and f = 4Hz (i) t = 3.1 iec, - 


1 (ii) t = 220.4sec, Rj *34. 
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Figure 8:The evolution of experiment B for D, = 2.0. The grid wss 
positioned 5 an below ihc interface and oacUlaied al f = 4H* 2 

cm atioke: (a) t = 5.2iec. R, - 5.0. D » 1.80 (b) i = g.l sec, R ; - 

M. D - 1.00 (c) l = 13.2 tec. R ; « 16.1. D - 0.78 (d) t = 91.9 »ec, 
Rj - 30.3, D - 0.54. 


Figure 5: Log plol of linear entrainment rue E = w e )w as a function 
of Richardson number Ri - gApJ/pei*. Grid loc ,u on , „ d 
frcquency: interface, f = 4Hi (o); 10cm above interface, f = 4Hz (,); 5 

an below interface, f = 4Hz (+) and f = 6Hz (x). For comparison, E A 
Hopfingers (1986) dau (solid line) with the slope of - 1.5 are 
included. 
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Figure 9:Buoyancy reversal parameter as a function of Richardson 
number: (a) Experiment A : +. Zo = 11 cm. f = 3Hz, S = 2 cm for 
Di = 0.7, 1.0 and 1.7; o, Zq = 10 cm, f = 4Hz, D; = 2.5. (b) 
Experiment B : o. Zo = 5 an. f = 4Hz. S = 2 an for Di =03 
0.8. 1.0 and 2.0. 
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Figure 6:The evolution of experiment A for Dj = 1,0 . The grid waj 
positioned 11 cm above the initial interface and oscillated at 3Hz with 
2 cm stroke: (a) t = 4.8»ec, Rj = 96.6; (b) t = 493.1 sec, R t = 38.5, 
D = 0.83; (c) t = 1053.4 *ec, Rj = 22.1, D = 0.59; (d) t = 1295.2 tec. 




Figure 10:Log plot of nonlinear entrainment rale E as a function of 
Richardson number and buoyancy reversal parameter D. Where Zq = 5 
cm, f = 4Hz and S = 2 cm for both runs with twe different initial values 
of D. Dash line m E v.i. Rj for linear case. 
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